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We analyze the redshift suffered by photons originating from an external source, traversing a 
collapsing dust cloud, and finally being received by an asymptotic observer. In addition, we study 
the shadow that the collapsing cloud casts on the sky of the asymptotic observer. We find that the 
resulting redshift and properties of the shadow depend crucially on whether the final outcome of 
the complete gravitational collapse is a black hole or a naked singularity. In the black hole case, 
the shadow is due to the high redshift acquired by the photons as they approach the event horizon, 
implying that their energy is gradually redshifted toward zero within a few crossing times associated 
with the event horizon radius. In contrast to this, a naked singularity not only absorbs photons 
originating from the source, but it also emits infinitely redshifted photons with and without angular 
momenta. This emission introduces an abrupt cutoff in the frequency shift of the photons detected 
in directions close to the radial one, and it is responsible for the shadow masking the source in the 
naked singularity case. Furthermore, even though the shadow forms and begins to grow immediately 
after the observer crosses the Cauchy horizon, it takes many more crossing times than in the black 
hole case for the source to be occulted from the observer’s eyes. 

We discuss possible implications of our results for testing the weak cosmic censorship hypothesis. 
Even though at late times the image of the source perceived by the observer looks the same in both 
cases, the dynamical formation of the shadow and the redshift images has distinct features and time 
scales in the black hole versus the naked singularity case. For stellar collapse, these time scales seem 
to be too short to be resolved with existing technology. However, our results may be relevant for 
the collapse of seeds leading to supermassive black holes. 


I. INTRODUCTION 

A widely studied collapse model in general relativity is the family of Tolman-Bondi spacetimes, which describe 
the collapse of a spherical dust cloud. The popularity of this model stems from the fact that in suitable coordinates 
the metric components and fluid fields can be expressed in terms of elementary functions, and thus provides exact 
solutions of the Einstein-Euler equations in the absence of pressure. Curiously, for a large class of initial data satisfying 
reasonable physical assumptions, these collapse models predict the formation of shell-focusing singularity a portion 
of which is null and visible to local observers, see for instance Refs. Furthermore, for suitable initial data in 

this class, part of this null singularity is visible even to observers in the asymptotic region; that is, the singularity 
is “globally naked”. While the Tolman-Bondi models, being spherical symmetric and having zero pressure, are 
highly restrictive from a physical point of view, still it should be emphasized that the occurrence of such globally 
naked singularities does not require any fine-tuning in the initial data within these models. Although there has 
been a significant amount of work regarding the properties of such singular spacetimes including the stability of 
the associated Cauchy horizon HIHHIH], relatively little attention has been placed on the manner that a globally 
naked singularity resulting from the complete gravitational collapse interacts with the rest of the Universe. However, 
nowadays, technological advances reach to the level where even the minutest predictions of general relativity can be 
placed under observational scrutiny, the current gravitational wave detectors UMB and the Event Horizon Telescope 
project [22] being prime examples that will test the theory in its strong field regime. These developments suggest 
that it is pertinent to start analyzing observational signatures associated with naked singularities. Even if the cosmic 
censorship conjecture would be proven rigorously some time in the future, ultimately its acceptance as a fundamental 
law of nature would require to be checked in the real Universe. 

In this spirit, in a recent paper |23| (paper I hereafter), we have proposed that the asymptotically measured 
frequency shift suffered by photons traversing a collapsing cloud can serve as a tool capable of differentiating between 
the formation of a globally naked singularity and the formation of an event horizon; i.e. this frequency shift can serve 
as an instrument to test the validity of the weak cosmic censorship hypothesis in terrestrial observatories. Evidence 
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supporting this proposal came from the analysis of radial photons from an external source passing through a collapsing 
spherical dust cloud. As shown in [24] , the total frequency shift of such photons relative to asymptotic static emitters 
and observers is always toward the red. Moreover, the analysis in I revealed that the photon’s energy measured by 
the asymptotic observer is gradually redshifted toward zero in the black hole case, whereas whenever a globally naked 
singularity forms the redshift of these radial photons remains finite and their frequency shift exhibits an abrupt cutoff 
once they have been detected by the observer who has crossed the Cauchy horizon (see Figs. 2 and 3 in I). 

The main goal of the present paper is to subject the proposal in I to further testing by allowing the incidence 
radiation to have nonvanishing angular momenta. Within the geometric optics approximation, this radiation consists 
of a collection of photons which are generated at past infinity, propagate freely through the collapsing cloud and 
are eventually detected by an observer in the asymptotic region. The central focus of this work is to examine the 
properties of the frequency shift of these photons. As will become clear in the course of this paper, the inclusion of 
angular momenta has interesting consequences at both the computational and the conceptual levels, and furthermore 
it adds a new flavor to the proposal in I. At the computational level, the evaluation of the redshift becomes a more 
difficult task since it requires integrating nonradial null geodesics inside the collapsing cloud. Nevertheless, as we 
show in this article, the nonradial photon’s frequency shift exhibits a qualitatively similar behavior than in the radial 
case, provided the direction is close enough to the radial one: first, the frequency shift is always toward the red. 
Second, this shift gradually converges to zero in the black hole case, whereas there is an abrupt cutoff in the naked 
singularity case. At the conceptual level, and to our pleasant surprise, we find that the nonvanishing photon angular 
momenta introduces certain “dark” directions at the center of the image of the source in the observer’s sky. These 
directions shape the optical appearance of the source, and imply that the collapsing cloud casts a growing shadow in 
the observer’s sky. This shadow resembles the well-studied black hole shadow. 

Bardeen [5^ describes the black hole shadow in the following way: “As seen by the distant observer, the black hole 
will appear as a ‘black hole’ in the middle of the larger bright source”. Cunningham and Bardeen |26j and more recent 
works |27H32] studied this shadow for a number of important illuminating sources. These works show that the black 
hole shadow is determined by tracing backwards null geodesics from the observer’s frame back to the illuminating 
source.^ Because the black hole region absorbs incoming radiation, certain directions in the observer’s sky will remain 
dark and the collection of all these dark directions determines the shadow that the black hole casts.Clearly, as long 
as the black hole is an equilibrium state, then its shadow remains steady provided the external illumination remains 
steady as well. For the case considered in this work, the changing background geometry adds a transient character to 
the perceived shadow. Nevertheless, the back-ray-tracing technique employed in the studies of black hole shadows is 
applicable to our problem as well. Through this technique, we compute the frequency shift suffered by photons and 
also study the characteristics of the shadow. 

Although the optical appearance and shadow associated with the formation of an event horizon have been studied a 
long time ago (see for instance Refs. jSSHHj)) no particular attention has been placed on the effects that the formation 
of a naked singularity has on the eyes of an asymptotic observer.'^ It is a priori not clear whether the formation of a 
naked singularity casts a shadow. If it does cast a shadow, one might wonder what its characteristic properties are. 
We delineate these issues by comparing the shadow that the collapsing cloud casts when the end state is a black hole 
to that alternative scenario, i.e. the case of a naked singularity. We find that in the former case the photosphere 
plays a crucial role as has been pointed out long ago by Ames and Thorne [5H]. Photons around the photosphere are 
responsible for the optical appearance, while photons grazing the horizon are received by the asymptotic observer as 
exponentially redshifted. In contrast, when a Cauchy horizon forms, the “visible” part of the singularity plays the 
dominant role in shaping up the shadow. The role that the singularity plays in the formation of the shadow is traced 
to a particular property of nakedly singular Tolman-Bondi spacetime. In such a spacetime, nonradial, future directed 
null geodesics escape from the central singularity and reach an asymptotic observer. This property has been proven 
for the case of marginally bounded Tolman-Bondi spacetime by Mena and Nolan Because of the important 

role that these null geodesics play in the formation of the shadow, we have recently verified their existence for the 
case of a nakedly self-similar Tolman-Bondi model |45] . and in the present work we show that Mena and Nolan’s 


^ Current day millimeter-wave very-long baseline interferometric arrays such as the Event Horizon Telescope m are on the verge of 
observing this shadow for Sagittarius A*, the supermassive black hole in the center of our galaxy. 

^ In this work we assume that the black hole or the naked singularity results from the collapse of a bounded system. If, for instance, a 
maximally extended Schwarzschild manifold is considered instead, the possibility that an observer in the vicinity of future null infinity 
could receive light originating from the white hole singularity ought to be taken into account and in this event the shadow would have 
different characteristics. We do not consider this possibility in this work based on the instability of white holes as demonstrated by 

Eardley m- 

^ We mention here that other ultracompact objects such as gravastars also generate shadows and in fact they may mask a black hole 
shadow. For recent work concerning this challenge, see, for instance, Refs. [31H37]. 

A notable exception constitutes the work in Ref. m- We shall discuss further their approach and results in the conclusions section. 
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conclusion remains valid for the large class of bounded nakedly singular Tolman-Bondi models. For the case of a 
Cauchy horizon formation, and once the observer lies to the future of this horizon, by a combination of analytical and 
numerical techniques we separate the past directed null geodesics emanating from the observer into those reaching 
the singularity and those reaching the asymptotic source. At first, we find that photons emanating from the central 
singularity and reaching the observer are infinitely redshifted and thus the singularity is invisible to the eyes of the 
asymptotic observer. In addition, we estimate the redshift suffered by the important class of photons that originate 
from the site of the asymptotic source, traverse the collapsing cloud and eventually reach the observer. We compute 
this redshift as a function of the viewing angles as perceived in the frame of the asymptotic observer and thus we 
extend the analysis in I where the considerations were restricted to radial photons. However, beyond the redshift 
analysis, the structure of the shadow as perceived by the asymptotic observer adds a new important component to the 
proposal put forward in I. This component, with potentially observational consequences, is related to the diverse time 
scales required by the shadow to develop in the eyes of the asymptotic observer. By an analysis of the snapshots of 
the image as perceived by the observer, we find that while for the case of a black hole formation, a few crossing times 
are sufficient for the shadow to occult the source from the eyes of the observer, this is not the case once a Cauchy 
horizon forms. It takes considerably many more crossing times so that the source can eventually be occulted from the 
observer’s eyes. Although for the case of stellar collapse the crossing time is extremely short, the situation is different 
for seeds leading to supermassive black holes. These issues are discussed in the conclusion section of the paper. 

The structure for the remainder of this paper is the following. In Sec.|llj because of the importance of null geodesics 
in this work, we begin with a brief derivation of the relevant equations of motion describing null geodesics (with and 
without angular momentum) on an arbitrary spherically symmetric background spacetime, and we derive the necessary 
formulas for the redshift of photons as measured by a certain class of preferred observers. Next, we specialize to the 
case of Tolman-Bondi spacetimes, describing the collapse of spherical dust clouds forming a locally or globally naked 
singularity. Based on methods from the theory of dynamical systems, we show the existence of infinitely many radial 
and nonradial null geodesics emanating from the singularity and analyze their local behavior. In particular, we prove 
that photons emitted from dust particles lying arbitrarily close to the central singularity suffer an infinite redshift. 
Next, in Sec. m we specify the initial data in our collapse model, and we provide a sufficient condition on the 
initial compactness ratio of the cloud under which the collapse leads to the formation of a naked singularity instead 
of a black hole. In Sec. IV we present a detailed discussion for our ray tracing method which ultimately allows us 
to determine the redshift image perceived by an asymptotic static observer. One of the key steps here consists in 
the determination of the critical angle of incidence a, which separates the null geodesics emanating from the naked 
singularity from those originating from the asymptotic source. Finally, in Sec.|V]we present the redshift images for a 
dust cloud forming either a black hole or a globally visible naked singularity, and compare the two cases. Conclusions 
are drawn in Sec. |VH and the article ends with the inclusion of two appendixes. In Appendixj^the proof of a technical 
result needed for establishing the existence of nonradial null geodesics emanating from the singularity is given. In 
Appendix]^ we discuss the aspects of the shadow of a naked singularity as seen by a comoving, free-falling observer 
in the particular case of self-similar Tolman-Bondi collapse, and we provide an analytic formula for the critical angle 


II. NULL GEODESICS IN SPHERICALLY SYMMETRIC COLLAPSING SPACETIMES 

We begin this section by reviewing basic properties of null geodesics propagating on an arbitrary spherically sym¬ 
metric spacetime. In particular, we point out the existence of a preferred vector field X, the Kodama vector field [46j . 
which reduces to the asymptotically unit timelike Killing vector field in the Schwarzschild case. We also show that 
the total redshift of photons measured by observers moving along integral curves of X assumes a particular simple 
form. Next, we specialize our considerations to the case of Tolman-Bondi dust collapse, and analyze the null geodesic 
flow in the vicinity of the naked singularity. Using standard results from the theory of dynamical systems, we show 
the existence of infinitely many radial and nonradial null geodesics emanating from the central singularity, a result 
that will play an important role in the considerations of Sec. |IV[ Finally, we prove that photons emitted from dust 
particles lying arbitrarily close to the central singularity suffer an infinite redshift when they travel along such null 
geodesics and are eventually detected by an asymptotic observer. 


A. General framework 

Consider an arbitrary spherically symmetric spacetime metric of the form 




( 1 ) 
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where $, 4' and r are functions of the local coordinates (t, R), and where d and ip are angular coordinates on the two- 
sphere. We are concerned with null geodesics propagating on the spacetime described by the metric given in Eq. Q. 
Because of rotational invariance, we can assume without loss of generality that the null geodesics are confined to the 
equatorial plane d = tt/2. The Lagrangian describing the motion is given by 


^ , n ^ cm (^\ _ 1 


—e 


2<S>{t,R) 


dt , 




dR\ , f dip 


( 2 ) 


where A is an affine parameter. Since ip is cyclic, we have conservation of total angular momentum, 

d£ _ 2 

and the Euler-Lagrange equations yield the first-order system 

dX 

^ = $e-2^(7r*)2-4/6-21' 
dX H 

dR _ 2-ifR 
dX 

^ = _$'e-2^(7r‘)2 + ^'e-2i'(^«)2 + ^/, 
dX 


( 3 ) 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


where we have introduced the momenta tt* := e^’^dt/dX and tt^ := e^'^dR/dX and where the dot and the prime denote 
partial differentiation with respect to t and i?, respectively. These equations should be supplemented with the nullity 
constraint C = 0. 

Eor the following, we introduce the vector field 


X := 



( 8 ) 


and the Misner-Sharp mass function |47j m defined by 

1- — — g{dr, dr) = + e~^'^r'^. (9) 

r 

Up to a sign, the vector field X is uniquely characterized by the properties of being orthogonal to the invariant 
two-spheres, of leaving the surfaces of constant areal radius invariant, i.e. X[r] = 0, and of being normalized such 
that 


g(X,X) = -(^l-^). (10) 

The existence of this natural vector field has been noticed long ago by Kodama [3S] in a different context. It follows 
that X is timelike outside the apparent horizon r > 2m. Associated with this vector field is the energy quantity 

£ := -g(X, p) = (Utt* -|- rn^) , (11) 

where here 

dt d dR d dip d 2 * t >9 f> d £ d 

P=XX5( + XXaR + dXSX; =' ' a +' " Sfl + J/fXX 

denotes the four-momentum of the particle. An additional useful equation is the change of areal radius r along any 
null geodesics, given by 


dr 

dA 


= re -I- r'e 


(13) 
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Combining Eqs. (§, ([TT] and ([T^ and using £ = 0 yields the equation 




dX 


r J 


(14) 


In the vacuum case, m is constant and X reduces to the time-translation Killing vector field. In this case, the energy 
quantity S is constant along the null geodesic flow and Eq. (14) expresses the familiar effective equation describing 
the radial motion for null geodesics. 

Finally, for later use, we compute the frequency shift of photons with respect to observers moving along integral 
curves of X. This frequency shift is given by 


^obs _ S(U, p)o6s Yj _ ^ 

“ g(U,p)e ’ “ V-g(X,X)’ 


which is valid as long as Xo^g and Xe are timelike, such that both the observer’s and the emitter’s four-velocity U is 
well defined. Using Eqs. (10 1 and 0 we obtain 


^obs 

j ^ 2m 

e 

^obs 

Ve \ 

^ 2 m 1 

r \ obs 

£e 


(15) 


In the Schwarzschild case, £ is constant along the null rays, and we recover the familiar expression describing the 
gravitational redshift. For a photon traveling from past to future null infinity through a collapsing spherical object 
of finite radius, Eq. (15) yields the remarkably simple formula 


OO 

^OO 


p+ 

Cl 


(16) 


where (£j“) is the energy £ of the photons at the moment they exit (enter) the collapsing object and and 
refer to the frequencies at future and past null infinity, respectively, as measured by asymptotic static observers. We 
will use this formula in the next section. 


B. Null geodesics in the Tolman-Bondi spacetime 


For the particular case of dust collapse, it is convenient to choose (t, R) as synchronous comoving coordinates [48] 
adapted to the flow of the dust particles. For such coordinates we have in Eq. ([^, <1> = 0, t = t is the proper time 
along the trajectories of the dust particles, and = 1 / 7 (t, R) = r'{T, R)/y/l + 2E(R). The areal radius r(r, R) 

of the spherical dust shell labeled by R at proper time t is determined by the equation 


1 

2 


m{R) 

r 


E{R), 


(17) 


with E{R) the energy of the dust shell R, see Ref. |7| for notation and further properties of the functions r{T,R), r, 
r' and 7 . Hereafter, we require the conditions (i)-(viii) in Ref. [7] to hold. These requirements imply regularity of 
the initial data and the absence of shell-crossing singularities (among other reasonable physical properties), and they 
lead to the formation of a locally naked shell-focusing singularity. 

For the following, we use the nullity constraint £ = 0 in order to replace with the single function 


7r« 

n := 7—, 

TT^ 


(18) 


in terms of which 


. Kl 1 11^1 n 

r vi - n2 ’ 1 r vi - n2 ’ 


(19) 


where we have taken the absolute value of ^ in order to make sure that tt'’’ > 0 such that the trajectory is future 
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directed.^ In terms of this new variable, the equations describing the null geodesics are: 


dr 

dX 

dR 

dA 

dn 

dX 


1^1 1 
r Vl - n2 ’ 

\i\ VI + 2 E n 
V ? VI - n2 ’ 


1^1 


^ 1-02 


-V1 + 2E + 

r 



and the change of the areal radius r along the null geodesic is described by 


dr 

dX 


— . ^ ^ (r + Vl + 2En) . 

r Vl -n2 V J 


( 20 ) 

( 21 ) 

( 22 ) 


(23) 


For the following analysis, we reparametrize the geodesics in terms of the coordinate R. This change of parametrization 
is legitimate as long as 11 7 ^ 0. Under this change of parametrization, the above equations take the form 


dr 

Ir 

dip 

Jr 

dr 

dR 

dU 

dR 


r' 1 

Vi + 2 £;n’ 

sign(£) r' Vl-n 2 

r V 1 + 2 U n 


1 + 


vT 




+ 2EII 

1 

yTTJE 




11 
r n 


(24) 

(25) 

(26) 
(27) 


Note that these equations are also valid in the case of radial null geodesics, for which 11 = ±1. The last two 
equations form a closed system of nonlinear ordinary differential equations for the variables (r, H), since the coefficients 
E,r,r,r',r' appearing on the right-hand side of Eqs (26) and (27) can be expressed in terms of R and r only. 


C. Radial and nonradial null geodesics emanating from the singularity 


In order to analyze the behavior of the solutions in the vicinity of the central singularity {t,R) = (0,0), it turns 
out to be convenient to replace the areal radius r with the new variable 


X := 


2m{R) 


and the coordinate radius R labeling the dust shells with s := 1/R. Note that y = 1 characterizes the location of 
the apparent horizon, while 0 < y < 1 defines the region lying to the past of it which includes the central singularity. 
Equations (26) and (27) can be rewritten as the following nonautonomous dynamical system for the variables (y, 11): 


^ f x \ _ x(Pt V mi - ( 

ds[ll) ~ - ( X 2 (i?,y,n) 


R—1/ s 


(28) 


with the functions Xi and X 2 given by 

^3 - / i?^^CR) 


Xi(i?,y,n) 

^ 2 (i?,y,n) 


X 

2 


R + -^yiy2EiR)^A 


X 


,R 


1 - 


' l + 2U(i7)y-2 y 

1 + 2E{R) n 


-3R-R 


2 ^{R) 

c{R) 


(i-n^) 

.a 


V A(5k31.fi 


3(fi)»/>' 




n 


fifim|i,fiU« 


^ Note that in the radial case 11 = ±1 and £ = 0. In this case, the equations describing the null geodesics are simply given by 
dr = n-^dR/y. 
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where the functions c{R), A{y,R) and K{y,R) are defined by 


c(i?) 

q{R) 

K{y,R) 

A{y,R) 

9{Q,y) 

fix) 


2m{R) 


I RE{R) 

1 

2^f^&h{q(R),y) - -^^^g(q(R),y), 
Rq{R) 2Rc{R)^^^^ 


A gWV R9'iR) 

^ 2c{R) y/l-q{R)2y2 q{R) 


y^ 


\/i - a/T^-^V 2 

fiqy) - fiq) 


- ndi^’V)’ 


5 \/1 — + arccos(a;), 


with y := ^rjR, see Ref. [7] for more details and properties of these functions. 

We are particularly interested in the behavior of the solutions of Eq. (281 in the vicinity of the central singularity 
(r, i?) = (0,0). Notice that in terms of the above notation, 


X = R\/ c{R) 



such that a finite, positive limit of x as i? —)■ 0 implies r/R —>■ 0. Therefore, any solution of Eq. (28) for which x 
converges to a finite, positive value when i? —>■ 0 corresponds to a null geodesics emanating or terminating at the 
central singularity. In order to understand the limit i? —>■ 0 (or s —>■ oo) we first analyze the autonomous system which 
is obtained by setting i? = 0 in the right-hand side of Eq. (28), that is. 


±fx 

ds \ n 


= X{0,X,I1) = Aox^ 


X n) ^ 

(i-n^)(3x-^) 


(29) 


where we have used the assumption E{0) = 0 which follows from condition (i) in Ref. [7], and where Aq := 
A(0, 0)/(2c(0)^/^) > 0. There are two critical points of this system with x > 0, namely 


(Xi) III) — (1,1), (X2>n2) 



The linearization about these points (keeping i? = 0 fixed) is 


77X(0,xi,ni) = Ao 



-DA'( 0 , X 2 , 112 ) — Aq 




In the first case, both eigenvalues are negative and thus the critical point (xi,ni) is an attractor for the system (29). 
In the second case, the eigenvalues are ±-\/2(2/3)^/^Ao, and hence the critical point (x 2 ,Il 2 ) is hyperbolic with 
associated one-dimensional stable and unstable local manifolds. A phase portrait for the system (29) is shown in 

Fig- [3 

After understanding the behavior of the autonomous system, we return to the nonautonomous system defined in 
Eq. (|2^, which we rewrite as 


( X 
ds \ n 


x(o,x,n) + i?r(R,x,n)|^^,/^ 


( 30 ) 


Y{R, X, n) := ^ [X(i?, X, n) - X( 0 , X, n)] , r> o. 


where 
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FIG. 1: The phase portrait for the autonomous dynamical system defined in Eq. (291. Note that we use the variable l/II instead 
of n in order to have a better view of the phase space in the region 0 < 11 < 1. The arrows show the direction o f the v ecto r 
field, the blue solid lines are the stable and unstable local manifolds associated with the critical point (x 2 ,112) = (•\/2/3, •\/2/3) 
and the red dashed line is an integral curve corresponding to a radial light ray. As the figure suggests, the basin of attraction 
for the critical point (xi. Hi) = (1,1) is the region that lies below the stable manifold. 


It is not difficult to verify that for small enough e > 0 and a suitable open neighborhood U of (x,n) = (1,1) 
in the map X : {0,e) x U — )■ x,n) >-)• X(i?, x,n) defined in Eq. ( p^ is continuous. Consequently, 

Y : (0, e) X [/ —)■ {R, x, 11) i—F {R, x^ n) defines a continuous map, and it is not difficult to prove that this map is 

bounded. Now we resort to the following general result whose proof we give in Appendix for completeness. 

Theorem 1 (Stability theorem, cf. Theorem 4.5 in Ref. [49]) Let U C be an open subset and let X ■. U ^ 

I—>■ X{x) be a continuously differentiable map and x* G U a point such that 


X{x*) = 0 


and such that all the eigenvalues of the linearization DX{x*) : K™ —>■ K™ at x* have a negative real part. Let 
Y : (0, oo) X f7 —)■ M™, (s,a;) i—)■ Y(s,x) be a bounded continuous map. 

Then, there exists si > 0 and an open subset V C U containing x* such that any (maximally extended) solution 
x(s) of the Cauchy problem 


/ 3? = ^( 2 ^) + s > So, 

\ x(so) = Xo 


with (sqjXo) G (sijOo) X V exists for all s > sq and satisfies 


(31) 


lim x{s) = X*. 

s—>-oo 


Remark: x{s) = x* for all s > 0 is not necessarily a solution of the Cauchy problem (311 with xq = x*, since we 
do not assume that Y{s,x*) = 0. 

According to Theorem]^ there exists an open neighborhood V of the point (x, H) = (1,1) and a value Rq > 0, such 
that any solution (x(s),II(s)) of the nonautonomous system (28) with (x(s), n(s))|g^^^^^ € V converges to the point 
(1,1) as s —>■ 00 . Physically, this implies that any null geodesic passing close enough to the point (i?, x,n) = (0,1,1) 
emanates from the central singularity. In particular, we have shown that there exists infinitely many null geodesics 
emanating from the central singularity, with or without angular momentum. As it will turn out, this property will be 


of central importance to the main results of this paper and its implications will become clear in Sec. IV 
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Clearly, the result described in the previous paragraph does not provide a full picture for the behavior of null 
geodesics in a vicinity of the naked singularity, since it leaves open the possibility for the existence of null geodesics 
emanating from the singularity with asymptotic values for (x,If) different from (1,1). A full picture in the generic 
Tolman-Bondi collapse case analyzed here is beyond the scope of the paper. For a full qualitative analysis of the null 
geodesic flow in the particular case of a nakedly singular self-similar Tolman-Bondi spacetime we refer the reader to 
our recent work [45]; see also Appendix [b| 


D. Infinite redshift 

After having showed the existence of infinitely many null geodesics, with and without angular momentum, emanating 
from the central singularity, we analyze the frequency shift the photons undergo along these geodesics as measured by 
free-falling observers. In particular, we show that photons emitted from a region very close to the central singularity 
are infinitely redshifted. 

The frequency shift along outgoing radial null geodesics is given by m 


^obs 


= exp 


Robs 

J ^iT+{R),R)dR 

Ra 


(32) 


where R i—>■ (t+ (R) , R) parametrizes the outgoing radial light ray, and R = Re and R = Robs denote the location of 
the free-falling emitter and observer, respectively. It follows from the explicit expressions for 7/7^ derived in Lemma 3 
of [T7] that 7/7^ ~ —1/R? for small i?, along radial light rays emanating from the central singularity with x —t 1 as 
i? —>■ 0, and consequently, there is an infinite redshift as Re —t 0, in agreement with Ref. |2]. 

Regarding the frequency shift along nonradial null geodesics emanating from the central singularity such that 
(x, n) —>• (1,1) as i? —)■ 0, we have 


^obs 

I'p. 


g(u, p) 


obs 


g(u,p) 



(33) 


where we have used Eqs. (12) and (19). Note that for a free-falling observer at i? = Robs > 0, we have Vobs > 0 and 
Ilgjg 7^ 1) since the null ray is nonradial. However, as we consider an emitter on a free-falling trajectory located closer 
and closer to the central singularity, re —>■ 0 and He —>■ 1, and we obtain again an infinite redshift. 


III. INITIAL CONDITIONS FOR THE COLLAPSING CLOUD 


In this section, we briefly discuss the family of initial data used in our analysis. The initial data for the collapsing 
dust cloud consist in the initial density profile p{R) and the initial profile for the radial velocity v{R) (for a precise 
definition of this velocity see |7]). For simplicity, we choose time-symmetric initial data for the dust collapse, implying 
that v{R) — 0 for all i? > 0. To be in accord with the simulations presented in I (Ref. |2^), the density profile is 
chosen in the following form: 


p{R) = pc>^ 



0<R<Ri, 
R> Ri, 


(34) 


where > 0 is the initial areal radius of the cloud, pc > 0 is the central density, and oq is a dimensionless parameter 
restricted to the interval 0 < oq < 1 which characterizes the flatness of the initial density p. The other parameters are 
given in terms of Oq by Oi = —3(6—5ao)/7, 02 = (8—5ao)/3, and 03 = —(9—5ao)/ll. The initial density profile defined 
in Eq. (34) is continuous and monotonously decreasing on the interval [0,oo), and provided 0 < uq < 12/5, it satisfies 
the assumptions (i)-(viii) listed in Sec. II of Ref. [7| such as smoothness, boundedness, non-negative mass density, 
absence of shell-crossing singularities, absence of initially trapped surfaces, with the exception of -differentiability 
at the surface of the cloud.® In particular, the profile given in Eq. (34) is smooth at the center where indeed p'(0) = 0 


® Nevertheless, the density profile is twice continuously differentiable everywhere. 
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and p"{Q) < 0. As shown by Christodoulou [5], these conditions lead to the formation of a shell-focusing singularity 
which is visible at least to local observers. 


be 


The initial compactness ratio C := 2mi/i?i of a collapsing cloud with initial density given by Eq. (34) turns out to 


C — -7ri?Jpc (1 — uq + Oi + 02 + 03 ) • 


(35) 


In our numerical simulations we fix ag = 1; nonetheless other values in the interval (0,12/5) lead to the same 
qualitative results. Below, we make the particular choice C = 16j77 =: Cg which leads to the formation of a black 
hole, and further consider initial data with C < Co/2 which guarantee the formation of a globally naked singularity; 
see paper I. 


IV. BACK-RAY TRACING THE PHOTONS FROM THE OBSERVER 

In the spacetime generated by the data described in the previous section, we consider a particular static observer at 
some distance robs ^ 2 mi from the collapsing cloud who registers photons that have traversed the collapsing cloud. 
As in I, we assume that the received radiation is generated by external sources distributed uniformly in the asymptotic 
region. This radiation does not interact with any intervening matter nor with the collapsing matter, and for simplicity 
we ignore any radiation generated within the collapsing cloud. However, unlike in I, here we allow this radiation to 
have nonvanishing angular momentum. Thus, within the geometric optics approximation, photons composing this 
radiation are moving on null geodesics originating from the asymptotic source, going through the collapsing cloud and 
finally received by the detectors carried by our asymptotic observer. In the spirit of the proposal in I, our focus in the 
next sections is to estimate the total frequency shift V^jv^ of the photons perceived by the observer, where hereafter 
and stand for the frequency of emission and detection, respectively. In order to compute the frequency shift 
of these photons, we trace back their path along past-directed null geodesics intersecting the cloud. Each 
null geodesic is characterized by its impact parameter b and the proper time Tobs at which the observer measures 
its corresponding redshift. In Fig. we illustrate a possible nonradial photon trajectory ( 7 ) in the case where the 
collapse forms a naked singularity, and we also show the particular radial light ray ( 70 ) penetrating the cloud at the 
moment of time symmetry, which we use as a reference to define Tabs = 0 . 



FIG. 2: Conformal diagram for the spherical collapse of a dust cloud forming a naked singularity. We illustrate the world line 
of the distant observer and the spacetime trajectories of two light rays 70 and 7 . The particular radial light ray 70 is the one 
originating from the distant source which penetrates the cloud at the moment of time symmetry and reaches the observer at 
time Tabs = 0 . 7 is an arbitrary trajectory of a null geodesic with angular momentum. 

We divide the process of tracing back the null geodesics in three steps. The first step consists in tracing back the 
null rays from the observer to the surface of the cloud, and thus determining the areal radius rf of the surface of the 
cloud at the moment the ray penetrates it and its angle of incidence a (defined below). No numerical ray tracing is 
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required in this step since we are dealing with null geodesics in a Schwarzschild spacetime, whose equations of motion 
are integrable. 

The second step involves tracing back the light ray numerically from the surface of the cloud through the cloud itself. 
In the case where the object collapses to a black hole, the light ray reaches the surface of the cloud again in the past. 
However, in the naked singularity case, there is an alternative possibility that may occur, once the observer passes the 
Cauchy horizon. As we have shown in the previous section, there exist future directed null geodesics emanating from 
the naked singularity. Therefore, it is possible that the ray-traced null geodesic hits the naked singularity instead 
of extending to the distant source. A consequence of this property is that once the asymptotic observer passes the 
Cauchy horizon, he or she registers a smaller number of photons originating from the source. In this step, we determine 
numerically the critical angle a which distinguishes between null geodesics emanating from the naked singularity from 
those that arrive from the external source at past null infinity. 

Finally, in a third step, we vary the parameters Tabs and b characterizing the null geodesic. For each value of these 
parameters, we determine whether the null geodesic originates from the source or from the naked singularity. In the 
first case, we compute the total redshift the photons undergo, while in the latter case, as we showed in Sec. [HDl the 
photons are infinitely redshifted. 

In the following, we describe the details involved in realizing these three steps. 


A. Step 1: Tracing back the null rays from the observer to the surface of the cloud 


In the exterior region the null geodesics are described by Eq. (14) with co nsta nt energy £. Introducing the impact 
parameter b := ijS and rescaling the affine parameter such that 1, Eq. (14) assumes the form 


dr\^ If 2mi 


1 

52’ 


(36) 


where mi is the total mass of the cloud. This is the equation for a one-dimensional mechanical particle in the potential 
(1 — 2mi/r)/r^, whose maximum l/{27m\) is located at r = 3mi (the photosphere). In the following, we restrict 
ourselves to absolute values for the impact parameter | 6 | < 6 , lying below the critical value 


b = V^TOi, 

corresponding to the maximum of the potential. Therefore, when | 6 | < 6 , 1/6^ is larger than the maximum of the 
effective potential, implying that the null geodesic must intersect the cloud in its past, since otherwise it must have 
emanated from r = 0 in the Schwarzschild region which is not possible for an observer lying to the future of the 
reference curve 70 , see Fig.[^ 

Let r = rf denote the areal radius of the surface of the cloud at the event x where the null geodesic intersects it. 
We define the angle a of incidence at x in the following way: denote by 

d d Id Id 

dr dR r dv r sm if d(p 

an orthonormal frame at x adapted to a free-falling observer comoving with a dust particle at the surface of the cloud 
and chosen so that ei is parallel to the outward pointing radial vector field at the event x. Then, we define the 
incidence angle a by 


= A (eg — cos a ei -I- sin a 63 ), 


(37) 


where A > 0 is a positive constant and Px denotes the four-momentum of a photon which is restricted to move on 
the equatorial plane i) = -k Geometrically, a determines the angle that the orthogonal projection of onto the 
invariant plane spanned by ( 01 , 03 ) makes with the radial outgoing direction 01 . According to this parametrization, 
Px is tangent to the ingoing (outgoing) radial geodesics for the particular value a = 0 (a = tt). 


Comparing Eq. (37) with Eq. (12) we find (recall that for dust collapse $ = 0, e ''' = 7 and t = t) 


tt’’ = A, 77r^ = — A cos a, £ = Arsma, 


such that n = — cos a. Introducing these relations into Eq. (11) and dividing by i yields 


b = b{a, r) = 


r sma 


\/l -|- 2Ei -\- \ -|- 2Ei cos Oi 


(38) 
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where we have used the equation jr' = \/l + 2Ei, Ei := E{Ri), and the free-fall equation evaluated at the 
surface of the cloud R = Ri. Note that for the particular values a = 0 or a = tt this formula gives 6 = 0, i.e. the 
correct impact parameter for ingoing and outgoing radial light rays. Inverting the relation ( |38[ ) we obtain the two 
solutions 


br 


sma = 


cos a = 




2mi 


,2 + (,2 (Im + 2£,) 

^2 _|_ ^2 ( 2 ^ _1_ 2^i) V r ^ \ r J 


(39) 

(40) 


for a as a function of the impact parameter and the areal radius r = of the surface of the cloud at the event x. 
Evaluating Eqs. (391 and (40) at 6 = 0 yields sin a = 0, cos a = ±1 which shows that in this case it is the lower sign 


that determines the correct solution a = n. More generally, it follows from Eq. (13) that 


dr ^1-^(1-^) 
i + g(^+2Ei) 


/- /2mi / 6^ / 2mi A 

tVI + 2Ei-^ — + 2EiJl - ^ (^1-—j 


(41) 


Since ^/2rnl/V+~2]^ < + 2Ei and 1 — (6/r)^(l — 2mijr) < 1 it follows that the lower sign leads to an increase 

of r along A, and hence the corresponding light ray reaches the observer. Thus, the correct sign in Eqs. (39) and (40) 


is the lower one, and the incidence angle a is uniquely determined by these equations and lies in the interval' 
7r/2 < a < 37r/2. In particular, this implies that a segment of the null geodesic in the immediate past of x lies inside 
the cloud, as required. 

Now that we have managed to express the incidence angle a in terms of the impact parameter 6 and the radius 
of the cloud when the null ray exits it, we need to establish the relation between rf and the proper time Tots at 
which the ray reaches the observer. We fix the time origin using the same normalization as in paper I. We choose this 
origin such that Tots = 0 corresponds to the moment at which the observer encounters the radial light ray 70 that 
entered the cloud at the moment of time symmetry (see Fig. [^. In order to compute the relation between Tots and 
rf we work in standard Schwarzschild coordinates (t, r) and first note that 


V 2mi , . 

1 [t ^o) — i ^0? 

robs 

where tg is the Schwarzschild time at Tabs = 0, and where in the last step we have used the assumption robs S> 2mi. 


Next, we use Eqs. (11) and (36), from which 


t-tt = 


dr 


N[r)^l-f,N{r) 


where is the Schwarzschild time at which the light ray exits the cloud and where for notational simplicity we have 
set iV(r) := 1 — 2mijr. Next, we use the fact that rf and are related to each other via the free-fall equations 

N{r)i = £1 = VI + 2E1, 


1 -2 "T-l 722 

o’’-= 

2 r 


where we recall that the dot denotes differentiation with respect to the proper time of a free-falling observer comoving 
with the dust particles. From these two equations we find, upon integration from the moment of time symmetry, 
where t = 0 and ri = i?i, to the point where the cloud has collapsed to radius rj*", that 


4 = + 2Ei f -^ 

4 N{r)^2E, 


+ 


2mi 


^ Notice that by restricting our considerations to the open interval 7r/2 < a < 37r/2 rather to 7r/2 < a < we have eliminated 

radiation that is incoming from infinity, grazes the surface of the collapsing cloud and returns back to infinity. For a discussion and an 
alternative description of this “backward emitted radiation” see Ref. m- 
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Combining the above equations, we obtain the following expression 


|&|,r+) = J -+ v/r + ^ f - ^ 

i N{r)Jl-^N{r) i N{r)j2E^ 

T'x y V 


2 mi 


(42) 


for the Schwarzschild time t at the moment the null geodesic intersects the observer’s world line, assuming that the 
geodesic has impact parameter b and exited the cloud at radius rf. By definition, 


to = t{Ei,0,ro), 

with rg the radius of the cloud at the moment the particular light ray 70 exits it. Hence, introducing the function 

1 


Fb(r) := 




- 1 , 


we obtain 


-<. = / 


2mi 


(43) 


The last two integrals on the right-hand side can be computed analytically, see Appendix C in [7]. The first integral 
converges as robs —t 00 , since T&(r) decays as 1/r^. Finally, we note that we can rewrite 


Pb{r) = 


N{r) 




so that the first integral converges as —>■ 2 mi. 

Gathering the results and taking robs 00 we finally obtain 

T{Ei,\blr+) 


where 

U(v) = i (\/«f - I>ly’ - “P {li - 

with oi := bi := ^/—2Ei and y := \JrjR\. In order to compute the remaining integral, it is convenient 

to introduce the new variable s := \b\/r in terms of which 


2bl 


arctan 




dr 


^l-^fV(r) (l + ^l-bA^N{r)\ 


4m 1 log 


U(yt)J 


(44) 


b^ 


dr 


\b\/r+ 


= \b\ 


ds 


sJl-'^Nir) (^1 -h yi-^fV(r)^ { 




(1+.yrr 


2 mi 

161 


The results of this subsection can be summarized as follows: if 7 is a null geodesic having an impact parameter b 
restricted to |5| < b, intersecting the observer’s world line at time Tabs > tq, then this geodesic 7 must intersect the 
cloud. The areal radius rj*" and the angle of incidence a [defined in Eq. (37l] at the event x at which the null geodesic 
exits the cloud can be determined from Tobs and b by first inverting the equation t{Ei, | 6 |,r[*") = Tobs to find rf and 


then using Eqs. (391 and (40) with r = r^ to determine a. Given r^ and a we can continue tracing back the null 


geodesic into the cloud, and this procedure is described next. 













































14 


B. Step 2: Tracing back the null rays inside the cloud and determining the critical angle a 


Starting from the surface of the cloud at the event x with r = rj*", at first we introduce the special angle a := a{b,r^) 
obtained from Eqs. (39) and (40) and subject to the constraint a € (7r/2,7r]. As explained in the previous subsection, 
the light rays with incidence angle a € (a,27r — a) are those that have large enough energy to make it over the 
potential barrier, see Eq. (36). Subsequently, for any a € (a,27r — a), we solve numerically the dynamical system in 


Eq. (28) in order to trace the nonradial null geodesic backwards into the dust cloud and determine whether it exits 


the cloud in the past or hits the singularity. If the collapse results into a black hole, only the former possibility can 
occur, whereas in the case of a naked singularity formation, the latter possibility may occur as soon as the event x 
lies to the future of the Cauchy horizon and the impact parameter b is sufficiently small in magnitude. As we have 
already discussed in the Introduction and have shown in detail in Sec. |IIC[ there exist future directed null geodesics 
emanating from the naked singularity, so it is possible that one of these geodesics intersects the surface of the cloud 
at the event x with incidence angle a within the admissible range {a, 2 tt — a) of the observer’s field of view. To deal 
with this new possibility, we introduce a critical angle a defined in such a way that null geodesics traced back from 
X hit the singularity in their past if a G (d, 27r — d) and exit the cloud in their past otherwise. If x lies to the past 
of the Cauchy horizon (in particular, this is always true in the black hole case), then d = tt and the light ray always 
exits the cloud in its past. However, as we will show below, as soon as the event x at the cloud’s surface crosses the 
Cauchy horizon, d decreases from tt to a value smaller than d as the radius of the cloud decreases to 2 toi, implying 
that there are more and more light rays in the observer’s field of view that emanate from the singularity instead of the 
distant source. Geometrically, the angle d determines the optical size of the naked singularity as seen by a free-falling 
observer which is comoving with the surface of the cloud. An analytic expression for d in the particular case of a 
self-similar Tolman-Bondi cloud is derived in Appendix [B} 

In the generic, bounded Tolman-Bondi collapse, an analytic determination of the critical angle d is probably not 
feasible. Therefore, we determine d numerically by bisection in a for every event x at the surface of the star such 
that Tobs > 0 and > 2mi. Figure]^ shows d and a as functions of 2mi/r^ for three different initial compactness 


ratios C (see the discussion at the end of Sec. Ill for the precise definition of C). In the left panel, C = Co, and 


thus a black hole appears, whereas the central and right panels correspond to the cases of naked singularities arising 
from collapsing clouds with C = Co/2 and C = Co/4, respectively. For every given r = Tobs and a discrete set of 
equal ly s paced val ues of the impact parameter b G [0,&), we numerically invert Eq. (44) to find r/" and then, from 
Eqs. (39) and (40) with r = r/", we obtain discrete sets of a S (d, d) producing the thin solid lines in each panel of 
Fig. ^blotice tHat the dependency of d on 2mi/r/" shows a kink at 2mi/r^ =2/3, corresponding to the moment 
the cloud’s surface coincides with the photosphere. The origin of this kink is related to the fact that the function 
1 — 6^(1 — 2TOi/r)/r^ > 0, appearing on the right-hand side of Eqs. (39) and (40), has a second-order root at r = 3 toi, 
implying a change of sign in the derivative of its square root. 





FIG. 3: The critical angle d (thick solid lines) as well as d (thick dashed lines) as functions of 2mi/r/. The thin solid lines 
correspond to slices of constant proper time of the observer. Left panel: Black hole case with initial compactness ratio C = Co- 
Central panel: Naked singularity with C = Co/2. Right panel: Naked singularity with C = Co/4, rch denotes the areal radius 
of the surface of the star at the moment it crosses the Cauchy horizon. Note that in the black hole case d = tt since the 
asymptotic observer always lies to the past of the Cauchy horizon in this case. 


In order to demonstrate the accuracy of our numerical method, in Fig. we present a self-convergence test for 
the numerical solution of the dynamical system in Eq. (28) computed by a fourth-order Runge-Kutta method with 
step size As = 0.002, for a null geodesic traced back from the event x until it intersects again the cloud’s surface 
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in the past. In Fig. [^we show results for an incidence angle a lying very close to the critical value d such that the 
corresponding null geodesic grazes the singularity. The plot shows that the relative numerical error decreases (left 
panel) by a factor ~ 16 (right panel) when the step size decreases by a factor 2 . 




FIG. 4: Self-convergence test for the numerical solution of the dynamical system in Eq. (28 1 for a null geodesic traced back 
from the event x at the surface of the cloud, passing very close to the singularity, and intersecting again the surface of the cloud 
in the past. The relative numerical error decreases (left panel) by a factor of ~ 16 (right panel) when the step size decreases 
by a factor of 2 . 


C. Step 3: Measuring the redshift of the photons originating from the source 


Once we have a light ray that exits the collapsing star in the past, we know that it must extend all the way to the 
past asymptotic region, where the corresponding photon is emitted with frequency v^. In this step, we evaluate the 
total gravitational frequency shift measured by the distant observer in the futu re a symptotic region, where 

denotes the detected frequency. This frequency shift is given by the simple formula (16 1 which only requires the 


computation of the quantity £ along the portion of the light ray inside the star. 

We numerically evaluate the frequency shift perceived by the observer at different proper times Tots > 0 and for 
different impact parameters b G (0,6) of the light ray.® To this purpose, we first determine the angle of incidence 
a G (a,7r) for each chosen value of Tabs and impact parameter b. In the naked singularity case we check whether 
a < a OY a > a. In the black hole case or in the naked singularity case with a < d, the light ray extends all the 
way to past null infinity according to the definition of a, and the total redshift is computed using Eq. (16 1 . In the 


naked singularity case with a > a, the light ray originates from the naked singularity and there is an infinite redshift, 
making it practically invisible. 

For the same cases as the ones presented in Fig. in Fig. we show the total redshift v^lv^ as a function of 
the viewing angle C (see definition below) as seen by the asymptotic observer for certain fixed values of Tabs > 0 . 
We denote by tch (Tsmi) the proper time at which the observer starts detecting light rays that exit the cloud inside 
the Cauchy horizon (photosphere). The viewing angle C is defined® in a similar way as the angle of incidence a, see 
Eq. (37). To define we choose the orthonormal tetrad 


eo = 


d 


yjN{robs) ^ '' ^ r Slu'd d^p 

adapted to the static observer, and we parametrize the photon’s four-momentum as 


= 'jNirobs) 


d 


e2 = 


1 d 


ea = 


d 


p = B (eo + cosCei + sinC 63 ), 


(45) 


® Because of spherical symmetry it is sufficient to consider positive values of 6, corresponding to the interval ol E (7r/2, tt] for the incidence 
angle. 

^ Here, a viewing direction is a spatial direction in the observer’s rest frame that describes the spatial direction of a photon originating 
from the source that has traversed the cloud and has been detected by the observer. It is related but not identical to the celestial 
coordinates of the observer’s sky. 
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with B > 0 a, positive constant. Using Eqs. 0 and p4| we find 


sin^ 


b 

T^obs 



2toi 

'f’ obs 


the angle C = 0 corresponding to photons with zero impact parameter (radial light rays). Since robs S> this 

relation simplifies to 



^ obs 


In particular, the largest viewing angle allowed by our assumption (| 6 | < 6 ) is ^ := b/robs- 


(46) 



FIG. 5; Total redshift measured by the asymptotic observer at certain proper times as a function of the normalized viewing 
angle ("/C of the collapsing object. Left panel: Black hole case with initial compactness ratio C = Co, for which T3mi/2mi ~ 10.8. 
Central panel: Naked singularity with C = Co/2. In this case Tchl^mi « 18.4 and Tsrnil^mi ~ 41.1. Right panel: Naked 
singularity with C = Co/4, for which Tch/2m\ ~ 47.5 and Tzmxl2m\ ~ 121.9. 

A comparison between the different panels in Fig. [^reveals that the redshift perceived by the asymptotic observer 
depends crucially upon the collapse outcome. In particular, regarding the total redshift suffered by radial photons 
(corresponding to C = 0) as a function of proper time, it can be seen from Fig. that our previous results in paper I 
are recovered. Namely, for the case of an event horizon formation, the left panel in Fig. shows that in an amount 
of proper time of the order lOmi crossing times, the observer finds that the ratio smoothly approaches zero. 

In contrast to this, for the naked singularity case, the central and right panels in Fig. [^show a sharp cutoff of the 
signal from the radial direction, which can occur even at crossing times of the order lOOmi, depending on the initial 
compactness ratio of the collapsing object. For a clearer picture of this behavior see also Figs. 2 and 3 in paper I. 

Furthermore, Fig. also shows that the behavior of the total redshift of radial photons persists for nonradial ones. 
That is, for every fixed viewing angle ^ 7 ^ 0, the redshift function decays smoothly to zero in the case of an event 
horizon formation, whereas in the naked singularity case the decay is slower initially but then shows a sharp cutoff to 
zero after the observer crosses the Cauchy horizon. Nevertheless, we see that as ( approaches </, the redshift seems to 
be independent of the collapse time. This can be understood by noting that C = C implies that we are dealing with 
null geodesics originating at the photosphere, and thus the associated redshift is independent of the collapse (for a 
related discussion of this point see Ref. [35] and in particular Fig. 2 of that paper). 

In addition. Fig. shows that in the naked singularity case the redshift behavior of nonradial photons as a function 
of C for late Tabs is also remarkable in the following sense: once the observer crosses the Cauchy horizon (robs > TCh), 
every light ray suddenly experiences an infinite redshift at some viewing angle C > Oj giving rise to a dark disk which 
will be discussed in detail in the next section. 


V. THE SHADOW OF THE NAKED SINGULARITY COMPARED TO THE SHADOW OF THE 

BLACK HOLE 

In this section, we provide a description of the “optical” appearance of the collapsing cloud as perceived in the 
eyes of an asymptotic observer, based on the implications of the graphs in Figs. [^ and We first examine the 
behavior of the angles a and d as a function of 2mx/r^. First, the panels in Fig.^show that the angle a exhibits a 
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qualitative behavior which is independent of whether an event horizon or a Cauchy horizon forms. This property can 
be understood by noting that as long as the surface of the cloud is in free fall, its behavior is entirely determined by 
the exterior Schwarzschild geometry. 

Next, we emphasize that in the case of an event horizon formation (see left panel of Fig. [^, any a £ [see 

Eq. (37) for the definition of a] defines a future directed null geodesic exiting the cloud and reaching the asymptotic 
observer, while in the past direction it traverses the cloud and reaches the source. Even though the observer is always 
allowed to receive photons from every direction such that a £ (djTr], after the surface of the cloud crosses inward 
the photosphere at r = 3mi, the bright directions in the observer’s sky progressively reduce due to the high redshift 
experienced by radial and nearly radial photons (see left panel of Fig. [^. In contrast to this, when a Cauchy horizon 
forms (see the central and right panels of Figs. [^and[^, the bright directions in the observer’s sky are determined 
by those values of a lying in the interval (d, a). Clearly, there is a dramatic reduction in the bright directions in the 
observer’s sky and the interval {a, a) shrinks to a point as the cloud’s surface reaches the photosphere.^*^ 

In order to provide the reader with an image of the observer’s view as he or she moves to the future, in Figs. ilZl an di 
we show color maps of evenly selected curves from Fig.[^ where the colors indicate the redshift’s magnitude, the darker 
regions corresponding to higher redshifts. Figure [^corresponds to the case of black hole formation. In this case, we see 
that the inner zone of the object progressively becomes darker, indicating a high redshift for nearly radial light rays. 
When the observer starts detecting light rays that exit the cloud inside the photosphere (that is, for Tobs > and 
< 3mi), there is already a high redshift zone at the center, and this zone becomes fainter and grows continuously 
until it covers nearly the whole region under consideration C < ^ as Tabs becomes large and rf -£ 2mi. These features 
show the dynamical formation of the shadow of the black hole, the final image consisting of a black region surrounded 
by a bright light ring due to the strong lensing effect at the photosphere. 

When a naked singularity results from the dust collapse instead of a black hole, we observe a completely different 
behavior of the shadow. This is mainly due to the nontrivial dependency of the critical angle a on the compactness 
ratio 2mi/ri; see Fig. As expected, a becomes strictly smaller than tt from the moment the observer starts 
detecting light rays that exit the cloud inside the Cauchy horizon (which corresponds to Tabs > Tch and rf < rch)- 
Moreover, the region a > a increases monotonously as the cloud collapses, and as a consequence, for r^bs > tch, the 
observer sees a dark disk concentric with the image of the cloud which occults the light rays coming from the source, 
see Figs. and This disk grows continuously until it covers the whole region under consideration at finite time 
Tabs = Tsmi (corresponding to = 3mi) which remarkably coincides with the intersection of a and d; see the central 
and right panels of Fig. This effect is also evident in the central and right panels of Fig. since there is a cutoff 
in the redshift function at some C > 0 for each Tabs > tch-, similar to the effect reported before in the case of radial 
light rays in paper I. Indeed, the increasing opaque disk arising in the naked cases corresponds to null rays coming 
precisely from the naked singularity, which are infinitely redshifted as discussed in Sec. |IID[ In this sense, we call this 
disk the shadow of the naked singularity. 



FIG. 6: Redshift snapshots at different proper times of the asymptotic observer as the dust star of initial compactness ratio 
Co collapses to a black hole. Here, rs^i/Smi « 10.8. 

In order to compare the dynamical formation of the different shadows discussed so far, in Fig. we plot their 
angular extension as a function of the proper time of the asymptotic observer. The dotted line corresponds to the 
shadow of the black hole, while the dashed lines correspond to the shadows of the naked singularities. An evident 


It may be noticed that the central and right panels of Fig. [^ show that after the cloud’s surface crosses the photosphere a < a. This 
reversal in magnitude is irrelevant from the point of view of the asymptotic observer since the trajectories of null geodesics with impact 
parameter b larger than b lie entirely in the Schwarzschild region a soon as rf < 3mi. 
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FIG. 7: Redshift snapshots at different proper times of the asymptotic observer as the dust star of initial compactness ratio 
Co/2 collapses to a naked singularity. In this case, Tc/i/2mi « 18.4, and T3mi/2mi « 41.1. The shadow of the naked singularity 
appears at Tabs = Tch, and then it covers the whole region at Tabs = Tsmi- 


nnnn 


oo 


T,J2m, = 1.1 T,J2m, = 15.0 T,J2m, = 30.0 T„,,/2m, = 45.0 T*/2m, = 60.0 r,J2m, = 75.0 



T*/2m, = 90.0 T„,,/2m, > T,J2m, 


FIG. 8: Redshift snapshots at different proper times of the asymptotic observer as the dust star of initial compactness ratio 
Co/4 collapses to a naked singularity. Here, rch/2mi « 47.5, and T3mi/2mi « 121.9. The same qualitative features as in Fig.j^ 
are present. 


difference is the time scale in which the shadows grow. It turns out that 97% of the shadow of the black hole develops 
in an interval of ~ 16.2mi crossing times, whereas the same growth of the shadow of the naked singularity with initial 
compactness ratio C = Coj2 takes ^ 45.4 toi crossing times, and ~ 148.8 toi crossing times for the naked singularity 
with the smallest initial compactness ratio C = Coj4. This means that the shadow of these naked singularities takes 
around 2.8 and 9.1 times longer to reach the same angular extension of the shadow of the black hole. 

Regarding the final state of the collapsing object and its appearance with respect to asymptotic observers, we note 
that a bright ring remains at late times in the black hole case, whereas the shadow of the naked singularity covers the 
whole region in finite time . Although the nature of the shadow in the naked singularity case is essentially different 
from the black hole case, since it is due to infinitely redshifted photons emanating from the singularity instead of 
highly redshifted photons originating from the illuminating source, it is important to stress that there is no difference 
in the image of the source for light rays with impact parameters larger than b at late times. This is due to the fact 
that at late times, light rays with b > b are scattered at the potential barrier, and thus they propagate entirely in 
the Schwarzschild exterior region. This suggests that observing the final state of the source’s image is insufficient to 
distinguish a naked singularity from a black hole. In this sense, the singularity remains censored by its own shadow 
and moreover, it could mimic a black hole shadow. However, as we have discussed, the dynamical behavior of the 
redshift during the formation of the shadow represents a possibility to distinguish between the two cases. 


VI. CONCLUSIONS 

This work has focused on an analysis of the redshift of a collection of photons illuminating a collapsing cloud and 
received by an asymptotic observer. The incident radiation is generated with arbitrary angular momenta by sources 
far away from the collapsing cloud, which for convenience are taken to be uniformly distributed. Moreover, this 
radiation does not interact with matter, but propagates freely, and eventually, part of it is detected by an asymptotic 
observer. Our results confirm the hypothesis put forward in paper I (Ref. [23)1. namely that the redshift can be 
considered as a tool capable of differentiating whether the outcome of a complete gravitational collapse is a naked 
singularity or a black hole. This confirmation consists in recovering the results in I for the case of radial photons, 





19 





FIG. 9: Normalized angular extension of the shadow of the black hole (dotted line) compared to the shadow of two naked 
singularities (dashed lines), as functions of the observer’s proper time. In the black hole case, the shadow becomes apparent 
around Tabs = T 3 mi (rf = 3mi), and then it grows continuously and takes an infinite time to cover the object, leaving a thin 
bright ring around it. In the naked singularity cases, the shadow appears in Tots = tch, and then it grows continuously and 
covers the whole region under consideration in finite time Tots = Tzmi ■ 


but fortifies the proposal by demonstrating that the redshift of certain classes of nonradial photons exhibits similar 
behavior as the radial ones. Our results show that photons in the incident radiation whose directions (as specified 
relative to the asymptotic observer’s frame) are close to the radial direction exhibit the same behavior as the radial 
ones. Specifically, there is a smooth decay to zero of the frequency shift for the case of event horizon formation, 
whereas in the naked singularity case the decay of the frequency shift is slower initially but shows an abrupt cutoff 
to zero after the observer crosses the Cauchy horizon. 

Beyond these new insights offered by the redshift of nonradial photons, the inclusion of angular momenta added a 
new element in the proposal in I. This new element consists in the structure and the properties of the shadow that 
the collapsing cloud casts into the sky of an asymptotic observer. Although the final state of the shadow turns out to 
be the same in both cases, we have found that the time scale required for the shadow to occult the source from the 
eyes of the asymptotic observer depends sensitively on whether the outcome of the collapse is a black hole or a naked 
singularity. 

Summing up, in this work we have presented new evidence to support the view in I that the structure of the 
redshifted radiation combined with properties of the shadow can be a tool in checking the validity of the weak cosmic 
censorship conjecture. 

The proposal put forward in I and tested further in this work is not the only proposal that has been considered in 
the literature. An alternative proposal is based on strong gravitational leasing effects. For instance, in Refs. [SUl - fS^ 
the gravitational leasing of sources due to a static, spherically symmetric solution of the Einstein-massless scalar field 
equations are compared to those of a Schwarzschild black hole. However, the choice of the nakedly singular spacetimes 
made in these works may be questioned on the grounds that no physical mechanism is known for the formation of 
these singular spacetimes starting from a regular set of initial data. In Refs. [HI 132, the effects of a naked Kerr or 
Kerr-Newman ring singularity are compared to those of a Kerr or Kerr-Newman black hole. Although one might 
think that such singularities could be created by throwing highly spinning particles into a nearly extremal Kerr black 
hole, all the results so far suggest that this scenario does not work due to spin coupling effects [55H5S] or backreaction 
effects (see Ref. m and references therein). In that sense, the above mentioned alternative applications to test the 
weak cosmic censorship lack physical relevance. 

One of the positive aspects of the proposal in I is that the choice of the Tolman-Bondi model as a background 
spacetime avoids the above mentioned difficulties. Whenever the collapsing cloud is considered to be of finite spatial 
extent, an exterior Schwarzschild region joins smoothly across the cloud’s surface. Such spacetimes can be generated 
by a regular set of initial data (for a discussion and properties of such data, see [71 [58]), and the maximal future 
analytic extension of the initial data surface can be determined. Moreover, by choosing the initial data appropriately. 
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we obtain spacetimes with prescribed causal structure. This property combined with the fact that the metric in 
the maximal future analytic extension of the initial data surface is known, allows one to perform an analysis of the 
behavior of causal geodesics through the collapsing cloud. Two recent works [SS] and [13] also employ a Tolman- 
Bondi background to test the validity of the weak cosmic censorship or to study the shadow of the singularity. In 
the first work, the consideration of the authors were restricted to marginally bound collapse. Unlike our scenario, 
they considered a situation in which the photons are emitted from dust particles inside the cloud. In contrast to our 
results, they did not find any specific signature in the emitted radiation that can differentiate between the formation 
of a naked singularity and the formation of a black hole. The work in [13] treats only the globally naked self-similar 
Tolman-Bondi collapse and considers the impact of nonradial null geodesics escaping from the central singularity on 
the shadow as the central question. Although their approach has some resemblance to the present work, they have 
assumed that the photons originating from the central singularity have unbound frequencies so that they are received 
with nonzero frequency by the asymptotic observer. Their choice of infinite frequencies invokes plausibility arguments 
based on the unknown theory of quantum gravity, and their conclusions differ from those obtained in the present 
work. 

It is interesting to mention briefly that a few conclusions reported in this work may extend beyond the family 
of the Tolman-Bondi spacetimes. For instance, it is well known since the work of Ames and Thorne [33] that the 
shadow arising during the formation of a black hole depends crucially upon the existence of the photosphere of the 
Schwarzschild exterior and has a weak dependence upon the details of the collapse. However, as far as the shadow 
of a naked singularity is concerned, many issues need to be considered in detail. For instance, it is not clear whether 
the existence of future directed null geodesics emanating from the central singularity is a generic property of null 
singularities or if it is associated with a particular feature of the Tolman-Bondi metric. If the first possibility occurs, 
then maybe some of the detail characteristics of the shadow can be extended to more general collapse models provided 
we have regular initial data whose evolution leads to a globally nakedly singular development. It is needless to say 
that these issues need to be checked much more thoroughly. 

What are the prospects that the proposal in I in combination with the results of this work, can be called to decide 
between the two competing alternatives in a realistic collapse event? 

As it has become clear from the present work, the crossing time U ■= GM/c^ = 10~^{M/ Mq)s with Mq = 2 x lO^^g 
the solar mass is the relevant time scale. Since typical stellar masses vary in the range (1 — IOO)M 0 , the corresponding 
crossing time U is of the order of a few ( 10 “^ — 10 “®)s which seems to be too short to be resolved with existing 
technology. However, cosmological scenarios offer better testing grounds. It seems by now that there is a consensus 
that supermassive black holes with masses in the range between IO^Mq and IO^Mq are the objects that energize 
quasars and active galactic nuclei. Even though this hypothesis resolves one problem, the same hypothesis raises 
the delicate question of the mechanism (or mechanisms) generating their enormous masses (for an overview of this 
delicate problem and further references see |60j 1. A currently popular scenario regarding their formation is based 
on the old idea of supermassive stars proposed long ago by Hoyle and Fowler [3IJ |33] (for an introduction of these 
configurations see Ref. [63]). In this scenario, at the early Universe, i.e. at times corresponding to cosmological 
redshifts 5 < z < 20, the dark halos led to the formation of supermassive stars with masses between IO^Mq and 
10^Mq, and these supermassive stars collapse gravitationally to form a black hole. These seed black holes subsequently 
grew via the accretion process to reach the range (10® — 10^)Mq. Although many assumptions enter in this scenario, 
it is important to realized that the validity of the cosmic censor conjecture plays a central role. Suppose, for example, 
that the outcome of the gravitational collapse of a supermassive dark star is a naked singularity instead of a black hole. 
In this event, the resulting shadow offers the best case scenario to be resolved by terrestrial instruments. However, the 
fact that the event takes place at a high cosmological redshift adds technicalities regarding the observability of such 
a collapsing event at least within the current technological capabilities. On the other hand, the analysis presented in 
this work stands on a firm basis and at the theoretical level differentiates between the two possibilities. 
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Appendix A: Proof of Theorem [T] 

The proof is based on standard arguments from the theory of ordinary differential equations; see for instance 
Refs. [49l[6l]- Without loss of generality we assume that x* = 0. Since U is open there exists (5 > 0 such that the 
ball Bs{0) of radius 6 centered at the origin is entirely contained in U: Bs{0) C U. Let e G (0, 6), and let so > 0 and 
xo G Bi;{0). Let a;(s) denote the maximally extended solution of Eq. , and define 

Sg := sup{s > Sg : x{s') G Bs{0) for all Sg < s' < s}. 

We claim that Sg = +oo provided that S > e > 0 are chosen small enough and sg > 0 is large enough, and that in 
this case lim x(s) = 0. The statement of the theorem then follows by choosing si large enough and V := B^{Q) with 

s—>oo 

e > 0 small enough. 

To prove the claim, we decompose 


X{x) = Ax + g(x), 

with A := DX(0) and g(x) := X(x) — Ax. According to Duhamel’s principle, x{s) satisfies 


S 

n(s) = + J 


1 


g(x(s')) + -y(s',x(s')) 


ds'. 


(Al) 


for all Sg < s < Sg. Since all the eigenvalues of A have a negative real part, there exist positive constants C > 0 and 
/3 > 0 such that 

|e^‘| < Ce-^* 

for all t > 0. Furthermore, by the differentiability of X{x) we can choose 5 > 0 small enough such that 

\9iy)\ < 

\y\ - 2 C 

for all y G ^^(O) \ {0}. Further, since Y is bounded, there is a constant K > 0 such that 

\Y{s,y)\<^ 


for all s > Sg and all y G Bs{0). It then follows from Eq. (Al) that 


S 

|a:(s)| < C'e-^(*-^«)|xg| + j ds' 


for all Sg < s < Sg. Hence, the function (j)[s) := ®°)|x(s)| satisfies 


S 

(j){s) <Ce+^ J ^(/)(s') + ds' =: ?/'(s). 


Sg < s < Sg 


Since 


we obtain 


f (^(s) + ^ , 


_d 
ds - 


= ('5-SO 


V(s) 


~ 2 5 


Integrating and using '0(sg) = Ce and |a:(s)| = e ®‘’^(/)(s) < e finally yields 

S 

|a;(s)| < = (®-®')ds' 

2 J s' 


(A2) 
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for all sq < s < Sq. The right-hand side is bounded from above by the constant Ce -I- K/sq. Therefore, if we choose 
e S (0,(5) small enough and sq > 0 large enough such that Ce + K/sq < S, it follows that the solution x{s) G Bs{0) 
stays forever inside the ball Bs{Q) and thus exists for arbitrarily large times s > Sq. 

In order to conclude the proof, we need to sho w tha t |a:(s)| converges to zero when s —>■ oo. For this, we first note 
that the hrst term on the right-hand side of Eq. (A2) vanishes for s —>■ oo since it is exponentially damped. As for 
the second term, using L’Hopital’s rule we find 


lim — [ ^ds' = lim —■ — --g- -- = lim — 

s—>-oo ‘2 J s—>-c 

So 


2 ^ ■ 


e 2 ‘ 


= 0 . 


This concludes the proof of the theorem. 


Appendix B: Appearance of the naked singularity in the self-similar case 


In this appendix, we compute the critical angle of incidence a for the particular case of a free-falling comoving 
observer in a nakedly singular, self-similar Tolman-Bondi spacetime. In this case, the presence of the homothetic 
Killing vector held ^ associated with the self-similarity allows one to analytically integrate the null geodesic flow. For 
a complete analysis for the behavior of the null geodesics in these spacetimes, we refer the reader to our recent work 
in [45j . The analysis and notation used in this appendix are based on that work. For related numerical work we refer 
the reader to Ref. [43| . 

For a self-similar collapsing spacetime, the metric coefficients in Eq. 0 are given by 

= 1, = F{x), r = RS{x), x — (Bl) 

II A 

with the functions 

^ 

P{x) ■■= S'(a;) := (1 - a; < 1/A. 

Here, A > 0 is a positive parameter determining the compactness ratio 

2m{R) _ 4 2 
R ~ 9^ 

at time t = 0. Note that in this model, this ratio is independent of R. Since r' = F > 0, there are no shell-crossing 
singularities in this model; however, there is a shell-focusing singularity at x = 1/A. When A < A* ~ 0.638, this 
shell-focusing singularity is naked. The causal structure of the resulting spacetime has been studied a long time ago; 
see for instance Refs. [4^165117(1] . Recently [45], we provided a full qualitative picture for the null geodesic flow (with 
and without angular momentum). Our analysis makes use of the homothetic Killing vector field and the main 
features of this flow depend on the effective potential 


W{x) := 


2;2 _ F'^(^x) 
S^X) ' 


x<-, 


see Fig. 10 for a plot of this function for the parameter value A = 0.6. As can be seen from this plot (cf. Lemma 2 
in [45]), W has three simple zeros located at a; = J^, a = 0,1, 2, the values Ja satisfying Jq < 0 < Ji < J 2 < 1/A. 
Geometrically, these zeros correspond to the set of points where the homothetic Killing vector field is null. W is 
positive for x < Jq and Ji < x < J 2 and negative for Jq < x < Ji and J 2 < x < 1/A, and it has a local maximum at 
some point x = Xc G (Ji, J 2 ). 

As a consequence of the conservation of angular momentum £ and the presence of the conserved quantity C associated 
with the homothetic Killing vector field, the motion in the x direction is restricted to the set W{x) < with 

/3 := tjC an “impact parameter”. For the following, the critical impact parameter /3c := Ij\/W{xc) corresponding 
to the local maximum of W plays an important role. The main qualitative features of the null geodesic flow which 
will be used here can be summarized as follows (cf. Figs. 2 and 4 in [15]): 


(i) The surface x = Jq \s generated by ingoing radial null geodesics terminating at the central singularity (t, R) = 
(0, 0), while the surfaces x = Ji and x = J 2 are spanned by outgoing radial null geodesics emanating from the 
central singularity. The surface x = Ji describes the Cauchy horizon, and it is generated by the earliest radial 
outgoing null geodesics emanating from the singularity. 
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FIG. 10: A plot of the effective potential W{x) for the parameter value A = 0.6. 


(ii) For small impact parameters, |/3| < /3c, all the outgoing null geodesics registered by a free-falling observer at 
constant R originate from the far past when x < Ji and from the central singularity when x > Ji. 

(iii) For large impact parameters, |/3| > /3c, a free-falling observer at constant R “sees” the following: all the 
outgoing null geodesics registered by the observer emanate from the far past when x < xi{/3), no null geodesics 
are registered when xi{/3) < x < X 2 {P) and all the outgoing null geodesics registered by the observer emanate 
from the central singularity when x > 2 : 2 (/3). Here, Xi{f3) < X2(/3) refer to the positive roots of W{x) = l//3^. 

The impact parameter /3 is related to the angle of incidence a defined in Eq. (|37[) in the following way: 


F(x) ±xQr(x) , (3 x"^ — F'^lx) 

___ mn cv _ 

X F F{x)Qi 3 {x)' S{x) X F F{x)Qp{x)' 


(B2) 


where Qp{x) := a/1 — P‘^W{ x). In the radial case, /3 = 0 and thus cos a = =Fl and sin a = 0, showing that in this 
case the upper sign corresponds to outgoing null geodesics, such that a = tt. In order to determine the appearance of 
the naked singularity as seen by a free-falling observer at constant R, we have to find the critical angle a such that 
all the outgoing null geodesics registered by the observer with a € (d, 27r — a) emanate from the central singularity 
and all the other outgoing geodesics originate from the far past (see Sec. IVB). Based on the properties (i)-(iii) listed 
above we find the following result for d. 

First, when x < Ji there are no null geodesics emanating from the singularity, since x = Ji describes the Cauchy 
horizon. Consequently, d = tt for x < R. Next, consider a fixed value of x in the interval Ji < x < Xc- In this case, 
the impact parameter (3 is restricted by the ineq ualit y /3^ < 1/W{x). However, only those null geodesics which have 
/3^ < /3^ and correspond to the upper sign in Eq. (B2) emanate from the central singularity (see Eig. 4 in [IS]). Hence, 
the angles of incidence a registered by the observer corresponding to null geodesics originating from the singularity 
are those for which 


— 1 < cos a < — 


F{x) FxQp^x) 
X F F{x)Qi3^{x) '' 


Jl < X < Xc- 


Next, assume x lies in the interval Xc < x < J^- As in the previous case, the null geodesics with impact parameter 
/3^ < 0^ corresponding to the upper sign in Eq. (B2| emanate from the singularity. However, when /3^ > j30 the 
null geodesics corresponding to both signs in Eq. (B2) are also registered by the observer as outgoing, and they both 
emanate from the central sing ulari ty (see the right panel of Fig. 4 in [45]). Therefore, cos a varies over all those values 
on the right-hand side of Eq. (B2) corresponding to the upper sign for which /3^ < 1/W{x) [1 > Q^x) > 0] and all 
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those values corresponding to the lower sign for which < 13“^ < 1/W{x) [QpS^) > Qp{x) > 0]. Therefore, the 
angles a corresponding to null geodesics originating from the singularity are those for which 


— 1 < cos a < — 


F[x) - xQp^{x) 
X - F{x)QpXx) ’ 


Xc < X < J2- 


Finally, when x lies in the interval J 2 < x < 1/A, W{x) is negative, and hence /3 can take any real value. However, 
as in the previous case, only those null geodesics corresponding to the upper sign in Eq. (B2) or those with the lower 
sign and /3^ > /?/ emanate from the singularity. Summarizing, we thus have 



f tt, 


X < Ji, 


a{x) = \ 

1 arccos ( 

' F{x)+xQfi^(x}\ 
x+F(x)Qi3^ix) 1 

, Ji < X < Xc, 

(B3) 

1 

arccos ( 

F(x)-xQi3^{x) \ 

^ x-F{x)Qfi^{x) J 

, Xc <X < j. 



Notice that a is continuous at a; = Ji and x^- Moreover, it has a well-defined limit as x ^ J 2 , since 

F{x) - xQfj^x) ^ 1 
X^J2 X — F{x)Q3 (x) 1 I 02 


The critical angle d as a function of the compactness ratio along the trajectory of a free-falling observer at constant 
R is shown in Fig. for the parameter value A = 0.6. 



FIG. 11: The critical angle of incidence a versus the compactness ratio C := 2m jr along the world line of a free-falling observer 
at constant R. The interval shown corresponds to the interval < x < Jah- As in the previous figure, the parameter A is set 
to 0.6. We note that the behavior of a in the self-similar case resembles the behavior of a. in the bound collapse case, see the 
central and right panels of Fig. 


Although in the parametrization given in Eq. (Bl) the singularity already appears in the t = 0 slice, one can show 
that the spacetime has a regular center at i? = 0 for negative values of t. Furthermore, it is possible to match the 
self-similar spacetime described by Eq. (Bl I to a Schwarzschild vacuum exterior spacetime. This can be accomplished 
at a constant R hypersurface by satisfying the standard junction conditions; see Ref. EU and references therein. The 
apparent horizon in the resulting spacetime is determined by the equation 2m = r, and it coincides with the event 
horizon in the Schwarzschild exterior, while in the interior it is described by the constant x = Jah surface given by 
S{x) = 4A^/9, that is, 


Jah 
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Since F{Jah) — Jah = 4A^/9 > 0, it follows that Jah > J 2 and that x = Jah is a spacelike surface. For A = 0.6 we 
obtain the numerical values 

Jo ~ -1.027, Ji ~ 1.129, J 2 ~ 1.495, Jah ^ 1-560. 
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